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Weyl points and line nodes are three-dimensional linear point- and line-degeneracies between two bands. In 
contrast to Dirac points, which are their two-dimensional analogues, Weyl points are stable in the momentum 
space and the associated surface states are predicted to be topologically non-trivial. However, Weyl points 
are yet to be discovered in nature. Here, we report photonic crystals, based on the double-gyroid structures, 
exhibiting frequency-isolated Weyl points with intricate phase diagrams. The surface states associated with the 
non-zero Chern numbers are demonstrated. Line nodes are also found in similar geometries; the associated 
surface states are shown to be flat bands. Our results are readily experimentally realizable at both microwave 
and optical frequencies. 



Two-dimensional (2D) electrons and photons at the ener- 
gies and frequencies of Dirac points exhibit extraordinary 
features LSD - As the best example, almost all the remark- 
able properties of graphene are tied to the massless Dirac 
fermions at its Fermi level Topologically L8], Dirac 

cones are not only the critical points for 2D phase transi- 
tions but also the unique surface manifestation of a topolog- 
ically gapped 3D bulk. In a similar way, it is expected that 
if a material could be found that exhibits a 3D linear dis- 
persion relation, it would also display a wide range of inter- 
esting physics phenomena. The associated 3D linear point- 
degeneracies are called "Weyl points". In the past year, 
there have been a few studies of Weyl fermions in electron- 
ics . The associated Fermi-arc surface states, quantum- 
Hall-effect lll5ll . novel transport properties Il6ll and the re- 
alization of the Adler-Bell-Jackiw anomaly Ill7ll are also ex- 
pected. However, no observation of Weyl points has been re- 
ported. Here, we present a theoretical discovery and detailed 
numerical investigation of frequency-isolated Weyl points in 
perturbed double-gyroid(DG) photonic crystals(PhCs) along 
with their complete phase diagrams and their topologically- 
protected surface states. PhCs containing frequency-isolated 
linear line-degeneracies, known as "line nodes", and their flat- 
band surface states are also presented. Unlike the proposed 
Weyl points in electronic system thus far, our predictions in 
photonics are readily realizable in experiments. 

Before proceeding, we first point out one intriguing distinc- 
tion between the 2D Dirac points and the 3D Weyl points. 
2D Dirac cones are not robust; they are only protected by 
the product of time-reversal-symmetry(T) and parity(P, inver- 
sion). In 2D, Dirac cone effective Hamiltonian takes the form 
of H(k) = v x k x a x + v z k z a z ; this form is protected by PT (prod- 
uct of P and T) which requires H(k) to be real. Thus, one can 
open a gap in this dispersion relation upon introducing a per- 
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turbation proportional to o y 
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that is imaginary; 



for example, even an infinitesimal perturbation that breaks 
just P or just T will open a gap. In contrast, 3D Weyl points 
are topologically protected gapless dispersions robust against 
any perturbation. In 3D, Weyl point dispersions are governed 
by the Weyl Hamiltonian H(k) = v x k x a x + v y k y o y + v z k z a z , 
the Gy term can exist only when PT is broken: this is a nec- 



essary condition for the existence of Weyl points. Since all 
three Pauli matrices are used in the Hamiltonian, one has no 
possibility of constructing a term that can open a gap in this 
two-band-degeneracy of 3D periodic systems, thereby making 
a single Weyl point absolutely robust to perturbations. The 
only way to eliminate and create Weyl points is through pair- 
annihilations and pair-generations of Weyl points of opposite 
chiralities, which typically requires a strong perturbation. 

Our starting point for obtaining Weyl points is a PhC whose 
bandstructure exhibits a three-fold degeneracy (quadratic in all 
three directions) at the BZ center (F) in an otherwise com- 
plete bandgap. This pseudo-gap system was reported in a DG 
PhC 11811 . whose real space structure is shown in Fig. [T^. 
When the air-sphere is not introduced, its gapless bandstruc- 
ture is plotted in Fig. |3i[33]. As shown below, applying 
different perturbations to this structure will lead us to PhCs 
with frequency-isolated linear point and line degeneracies (i.e. 
Weyl points and line nodes). 

A DG PhC consists of two single gyroids (SG) 
in a body-centered-cubic(bcc) lattice. An SG surface 
is a bi-continuous triply-periodic minimal surface de- 
fined by g(r) = sin(27rx)cos(2;ry) + sin(27ry)cos(2;rz) + 
sin(2;rz)cos(27rx) lll9ll . The red gyroid in Fig. [T^ is defined 
by filling the inner space of the isosurface(g(r) > 1.1) to be 
high refractive-index material( y/e =n = 4) and air otherwise; 
magnetic permeability jl is unity everywhere. (This setting 
corresponds to Germanium and air at optical frequencies.) 
The blue gyroid is the inversion pair of the red gyroid with 
respect to the origin; the two gyroids do not overlap in space. 
The bandstructures of both SG PhC and DG PhC are plotted in 
Fig. [2^ by orange and blue colors respectively. The SG PhC 
has a 32% complete bandgap between the 2 nd and 3 nd bands 
from 0.42 to 0.58 in normalized frequencies fl8ll . The space 
group of the SG is 74i32 (No. 214) IMIHI]; the DG structure 
belongs to laid (No. 230), which is a direct-product group 
of 74i32 and inversion. Both of the space groups are non- 
symmorphic and have three-dimensional representations 1E2I1 . 
This explains the three-fold degeneracy in the DG PhC band- 
structure, among 3 rd , 4 th and 5 th bands at the center of the BZ 
(F), highlighted with green ellipses in Fig. [2^. The 1 st and 
2 nd bands are almost degenerate; so are the 3 rd and 4 th bands: 
they concave down and touch the 5 th band which concaves 
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wise complete frequency gap. It is worth pointing out that 
this bandstructure, although not exhibiting Weyl points, is in- 
teresting in itself in analogy to the line-node semimetals 112311 . 
We show its flat surface dispersions towards the end of this 
paper. 

In what follows, we break the PT symmetry to obtain Weyl 
points of photons for the first time. We start with the struc- 
ture from Fig. [2^. First, we break P while preserving T. Since 
T maps a Weyl point at k to — k with the same chirality, there 
must exist at least two other Weyl points, both of opposite chi- 
rality, to neutralize the whole system. So the minimal num- 
ber of Weyl points in this case has to be four. We break P 
by placing only one air s phe re on one of the gyroids (but not 
the other) at (|.-Q)fll34ll. as illustrated in Fig. Ola). Un- 
der this pure P-breaking perturbation! 35], two pairs of Weyl 
points, shown in Fig. [2k, emerge along T-N and T-H direc- 
tions. The fact that all the Weyl points appear along high- 
symmetry lines significantly simplifies the analysis. There are 
no other states in the vicinity of the Weyl points' frequencies. 

Second, DC magnetic fields (B), along different directions, 
are applied to the original DG PhC structure in Fig. [2^ to 
break the T while preserving P. We assume the high-index gy- 
roid material is gyroelectric and use a generic model ll24ll to 
describe its magnetic response. When B is along z, we assume 
the permittivity tensor takes the form of 
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FIG. 1 : Real-space unit cell and reciprocal-space BZ of the 3D 
DG PhC. a, Real space geometry in a bcc unit cell where ai = 

(-l,l,l)f, a 2 = (l,-l,l)f and a 3 = (1,1,-1)|. The two identi- w here det(e(|B|)) = (eft (|B|) - e} 2 (\B\))e = £ 3 ; this constant 
cal gyroid structures of red and blue colors are high refractive index 
(n = 4) materials; they are inversion pairs of each other with respect 



to the origin(o). An air sphere (r/a = 0.13) located at — 
breaks the inversion symmetry of the system, b, The BZ of the bcc 
lattice. Weyl points and linear line- degeneracies of investigation in 
this letter always lie in the green plane through the origin(r). T-N is 
along [101] and T-H is along [010](y). c, An air-isolated DG surface 
can be formed by terminating the perturbed gyroid (red) but not the 
other (blue). The SG PhC on the top has a large complete bandgap 
as shown in Fig. [2^. 



up. The three-fold degenerate point is well isolated in fre- 
quency from other parts of the dispersion diagram of the DG 
bandstructure, making it an ideal starting point for applying 
symmetry-breaking perturbations. 

The three-fold degeneracy of quadratic dispersions at T can 
be lifted by breaking the 74i32 space group without breaking 
P or T symmetries. This is done by replacing a part of the gy- 
roid material with two air- spheres (one on each gyroid). The 
first air-sphere is placed in the red gyroid at — ^, ^)a, as 
illustrated in Fig. [T^; and the other is its inversion pair in the 
blue gyroid (not illustrated in Fig. [T^). This perturbation lifts 
the 5 th band out of the three-fold degeneracy with the 3 rd and 
4 th bands at T, as shown in Fig. [2]b. The 4 th and 5 th bands 
linearly cross each other, forming a closed line-degeneracy 
around the T point in the T-N-P-H plane, inside an other- 



determinant condition ensures the dispersions as a whole do 
not move much in frequency with the external DC B fields. 
The dimensionless effective magnetic field intensity is de- 
fined as |B| = £12/ £ in this paper. When B field is along 
other directions, the corresponding £ tensor can be obtained 
through coordinate transformations. (Note the T-breaking can 
be equally well implemented via jl for gyromagnetic mate- 
rials [25].) Under this pure T-breaking perturbation, only a 
single pair of Weyl points emerges along the direction of the 
magnetic field. This is the minimum number of Weyl points 
that can exist under the inversion symmetry. These two Weyl 
points are still frequency-degenerate: P maps a Weyl point at 
k to — k with the opposite chirality. An example of this is 
shown in Fig. [2]i. 

Third, we apply both P and T breaking perturbations at 
the same time to observe the phase transitions between the 
two(II) Weyl points in the pure T-breaking phase and the 
four (IV) Weyl points in the pure P-breaking phase. Inter- 
estingly enough, different magnetic field directions produce 
strikingly different phase diagrams. When B is applied along 
T-H, only two phases exist: the T-breaking dominated phase 
(II) and the P-breaking dominated phase (IV). The pure P- 
breaking phase, shown in the contour plot Fig. [2t, has four 
Weyl points: two with positive chiralities along T-H and two 
with negative chiralities along T-N. Applying magnetic field 
along the T-H direction drives the two negative-chirality Weyl 
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FIG. 2: Gapless photonic bandstructures of the DG PhCs without (a) and with perturbations (b,c,d). a, The original DG PhC bandstruc- 
ture has a three-fold degeneracy at T among the 3 rd , 4 th and 5 th bands in a pseudo-gap. The SG PhC has a huge frequency gap covering the 
interesting frequency region of the DG PhC dispersion, b, We place the two air-spheres on the two gyroids: one air-sphere (r/a = 0.07) is 
located at — ^)a; the other is its inversion symmetric pair. Under this perturbation, the 4 th and 5 th bands linearly touch in a closed line 
around T in the T-N-P-H plane. The linear crossing line is highlighted by a green stripe; this structure does not exhibit any Weyl points yet. 
c, We apply a P-breaking perturbation by placing only one air sphere (r/a = 0.10) in one of the gyroids but not the other. Two pairs of Weyl 
points appear(highlighted by the green circles): one pair appears along T-H and the other pair appears along T-N. d, We apply a T-breaking 
perturbation (P conserves) by applying DC magnetic field (dimensionless |B| = 0.875) on the DG PhC without air- spheres. Only one pair of 
Weyl points appears (highlighed by the green circle) along the direction of the magnetic field (T-N). Since P and T are not broken at the same 
time in these PhCs, k and -k are degenerate in the bandstructures. A few lowest-value equi-frequency contours of the frequency difference 
between the 4th and 5th bands are shown (hexagonal insets) for each bandstructure in the T-N-P-H plane; the contour spacing is 0.004 in 
normalized frequency; "+" and "-" are used to label the chiralities of the Weyl points. All the dispersion behaviors close to the degeneracy 
points in this figure can be well described by the low-energy-theory model in the supplementary material. 



points towards the positive one along T-H. Increasing the 
magnetic field further annihilates two of the Weyl points of 
opposite chiralities and the system enters the T-breaking dom- 
inated phase which has only two Weyl points along T-H. An 
example where only T is broken is shown in Fig. [2]l Detailed 
description of this phase diagram is shown in the supplemen- 
tary Fig. A. 

An even richer phase diagram, shown in Fig. [51 appears 
when we switch the magnetic field to the T-N direction. The 
system undergoes two phase transitions from the P-breaking 
dominated phase(IV) to the T-breaking dominated phase(II). 
The extra intermediate phase(VI) of six Weyl points is gen- 
erated when one of the negative-chirality Weyl points in the 
pure P-breaking phase (red contour plot) splits (orange con- 
tour plot) into three Weyl points (blue contour plot) under the 
increasing B field. Among these three Weyl points, the origi- 
nal Weyl point flips its chirality from negative to positive; the 
other two negative Weyl points move away from their creation 
position towards the neighboring two positive ones. The two 
pairs of Weyl points eventually meet (green contour plot) and 
annihilate simultaneously, leaving the remaining two Weyl 
points (purple contour plot) along the direction of the mag- 
netic field(r-N). When two Weyl points annihilate, the band 
dispersions are quadratic along the directions in which they 
meet. 



In general, P-T phase diagrams depend on the form of the 
perturbations: the DC magnetic field can be oriented along an 
arbitrary direction; in addition, there are many different ways 
to break the inversion symmetry. 

Weyl points are topologically stable objects in the 3D BZ: 
they act as monopoles of Berry flux in momentum space, and 
hence are intimately related to the topological invariant known 
as the Chern number 12611 . The Chern number can be defined 
for a single bulk band or a set of bands, where the Chern num- 
bers of the individual bands are summed, on any closed 2D 
surface in the 3D BZ. The difference of the Chern numbers 
defined on two surfaces, of all bands below the Weyl point fre- 
quencies, equals the sum of the chiralities of the Weyl points 
enclosed in between the two surfaces In particular, we 
illustrate in Fig. [4^, the Chern numbers of the 2D planes per- 
pendicular to r-H(y) ; they are closed surfaces in the 3D pe- 
riodic BZ. The first contour plot (red), of a pure P-breaking 
phase in Fig. is used as an example. The Chern number of 
all lower bands on the plane is plotted in light blue on the left 
side. The Chern number vanishes when the plane is outside 
of the Weyl points. It changes by one when the plane moves 
across one Weyl point and changes by two when the plane 
crosses two Weyl points of the same chirality. 

As a result of the non-zero Chern number of all the lower 
bulk bands (1 st , 2 nd , 3 rd and 4 th ), there exists topologically- 
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FIG. 3: Phase diagram of the Weyl points when the magnetic field is applied along T-N. Each phase is labeled by the number(Roman 
numerals) of Weyl points. The phase transition line is found by scanning the B field amplitude in the calculation for three different P-breaking 
strengths. The pure P-breaking phase (x-axis) is highlighted in red; it is realizable at the optical frequencies with essentially lossless materials. 
Six 2D contour plots show the locations of the Weyl points at representative positions in the phase diagram; the contour spacing is 0.002 in 
normalized frequency. The contour plots share the same color as the circles of their corresponding positions in the phase diagram. Two surface 



dispersion plots show the entire band dispersions fi)(k) of the 3 rd , 
Fig.Eb. 



4 th and 5 th bands for two particular cases, where k lies in the green plane of 



protected gapless chiral surface states inside the bandgap (be- 
tween 4 th and 5 th bands) away from the Weyl points 112711 . In 
order to prevent the surface states from radiating into air, we 
interface the DG PhC with an SG PhC by removing only the 
perturbed gyroid above the (101) surface: this is illustrated in 
Fig. [it. The surface states of finite k y , are trapped by the par- 
tial gap of the DG PhC and the full 3D gap of the SG PhC. We 
define a termination parameter t(0 < t < 1) to indicate the 
periodically equivalent termination positions along the [101] 
direction, x = is set at the origin of the unit cell as shown in 
Fig. [Tfc. An example of the non-trivial gapless surface disper- 
sion, along a line cut in the 2D surface BZ in Fig. [4^, is plot- 
ted in Fig. @]b. A surface mode profile is also shown. Since 
T is not broken in this example, the surface dispersions are 
degenerate between k and — k, while Chern number flips sign 
under T symmetry. This means the surface modes have oppo- 
site chiralities when travelling along +y and — y directions at 
the same frequency. The above analysis on Chern numbers 
and surface modes applies for states along other directions 
and for other Weyl phases as well. We stress that the non- 
zero Chern numbers and the associated chiral surface states 
in this example are obtained without breaking T tl|,|25|]; thus, 
this phenomenon is readily realizable at optical frequencies. 

The surface states of the other two gapless PhCs in Fig. 
[2k and b do not involve non-zero Chern numbers, but they 



can exhibit very flat bands which is also an interesting fea- 
ture. Fig. [5fc shows one surface state of the DG PhC of 
quadratic point-degeneracy in the bulk pseudo-gap. When 
T varies periodically, the surface dispersion, at every surface 
k point, moves from the air-band(conduction band) through 
the pseudo-gap to the dielectric-band(valence band) lEsll . At 
the center of the BZ, the surface dispersion is pinned into the 
bulk states at the degeneracy point. The surface dispersion of 
T = 0.0, apart from the H-P region, is very flat. The high 
density of states associated with the flat surface dispersion 
can be potentially useful for enhancing the efficiency of light- 
emitters at the surface. Even more interesting surface states 
are shown in Fig. 0) for the line-node PhC. The nodal line 
bulk states project onto the surface BZ as a closed line which 
separates the surface BZ into two disconnected areas. So the 
surface dispersions can be flat bands in these two separate re- 
gions of the BZ. The green dispersion in Fig. Q) has all its 
frequencies nearly degenerate inside the line-node area; while 
the red dispersion is relatively flat in the rest of the BZ. The 
general features of the flat surface dispersions do not change 
when the line-node PhC is terminated by other means. For 
example, when the SG PhC is replaced by air, one could se- 
lectively enhance, by changing the surface termination, the 
light emission of surface sources into either radiative or non- 
radiative modes. 
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FIG. 4: Topological surface states of a Weyl-point PhC. a, Illus- 
tration of the Chern number of the lower bulk bands on the plane 
moving along k y in the 3D bulk BZ. b, An example of a topologi- 
cally protected surface state of a Weyl point PhC. The projected bulk 
states are plotted in blue in the surface bandstructure; the surface 
state dispersion relation is plotted in red. The spatial field distribu- 
tion (|E| 2 ) of one surface mode, at the center of the k axis, is plotted 
with the contour outlines of the refractive index profile. The surface 
state is localized at the interface between the DG PhC and the SG 
PhC. 
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FIG. 5: Flat surface dispersions of two pseudo-gapped PhCs 
without Weyl points, a, Surface dispersion at the interface between 
the SG and the original DG of a bulk quadratic point-degeneracy as 
in Fig. |2b. b, Two surface dispersions at the interface between the 
SG and the DG modified by two air spheres as in Fig. |2]b. They are 
flat in two different regions of the surface BZ separated by the linear 
line-degenerate bulk states. The surface BZ is colored in the region 
where the surface states are flat. 



In conclusion, this letter provides a detailed design of the 
minimum numbers of frequency-isolated Weyl points and line 




nodes along the high-symmetry lines and planes in the BZ of 
the DG PhCs. Complete P-T phase diagrams are investigated 
through annihilations and creations of the Weyl pairs. The 
topologically non-trivial surface states of Weyl-point PhCs 
and flat surface dispersions of line-node PhCs are presented. 
This work potentially opens doors to new paradigms in pho- 
tonics: topologically-protected chiral surface states realizable 
at optical frequencies for the first time, radiation-controllable 
photonic density of states enhanced by flat surface disper- 
sions, novel transmission properties of light in the vicinity 
of the degeneracy frequencies, and other unconventional phe- 
nomena associated with the density-of- states and dispersions 
of the close to the degeneracy points that were not available 
before. Experimentally, our design can be realized with many 
available techniques. Germanium can either be directly syn- 
thesized J29ll into DG strucutures or backfilled [30] into DG 
polymer templates made by block-copolymers Dill or inter- 
ference lithography J32I1 . The templates for geometrically- 
perturbed DGs can be made by 3D lithography or 3D printing. 
All presented results without magnetic fields are readily real- 
izable with Germanium at optical frequencies; Weyl phases 
involving T-breaking can be readily demonstrated with gyro- 
magnetic materials at microwave frequencies ll25ll . 
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FIG. A: Phase diagram of Weyl points when the magnetic field is applied along T-H. Each phase is labeled by the number(Roman numerals) 
of Weyl points. We start with the structure from Fig. 2a: the gray contour plot. The phase transition line is found by scanning the B field 
amplitude in the calculation for three different P-breaking strengths. The pure P-breaking phase (x-axis) is highlighted in red; it is realizable 
at the optical frequencies with essentially lossless materials. Six 2D contour plots show the locations of the Weyl points at representative 
positions in the phase diagram. The contour plots share the same color as the circles of their corresponding positions in the phase diagram. 
A few lowest- value equi-frequency contours of the frequency difference between the 4th and 5th bands are shown in each contour plot. The 
contour spacing is 0.002 in normalized frequency. "+" and "-" are used to label the chiralities of the Weyl points. Two surface plots show the 
entire band dispersions of the 3 rd , 4 th and 5 th bands for two particular cases of the contour plots. 



In Fig. [Al when B is along T-H, only two phases exist: the T-breaking dominated phase(II) and the P-breaking dominated 
phase(IV). The phase transition line is found by scanning |B| for different P-breaking strengths (radii of the air spheres). The 
annihilation and creation of Weyl points are shown by contour plots in a 2D BZ slice(green plane in Fig. lb), where all the 
Weyl points in this letter happen to be located. The pure P-breaking phase, third(from left) contour plot in red, has four Weyl 
points: two with positive chiralities along T-H and two with negative chiralities along T-N. Applying magnetic field along the 
T-H direction drives the two negative-chirality Weyl points towards the positive one along T-H (blue contour plot). On the phase 
transition line where the three Weyl points meet (green contour plot), the dispersions of the two bands are cubic(±|& 3 |) along 
the direction which they meet, but still linear in other directions in k space. Increasing the magnetic field further annihilates two 
of the Weyl points of opposite chiralities and the system enters the T-breaking dominated phase(II) which has only two Weyl 
points (purple contour plot). 
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MODEL OF LOW ENERGY THEORY 

In order to further elucidate the physics and to tie our work more closely to the Weyl fermions in electronics, we constructed 
a three-band effective Hamiltonian(//o) with three perturbation terms(Vi, V2 and V3) to model the 3 rd , 4 th and 5 th bands of the 
DG PhC system. This model reproduces qualitatively all the bandstructures in Fig. 2 and the phase diagram in Fig. El In order 
to reproduce the phase diagram of Fig. 3, a more intricate model is required. 

H(k) = flb(k) + Vi + V 2 (k) + v 3 
if (k) = a 1 |k| 2 -a 2 |k-L| 2 

y 1 = - 7l (|L.n p | 2 -l) (1) 
y 2 (k) = r2 (k-no)(L-n ) 

V 3 = j3(L-B) 

/o o\ /o -1 \ / 1 0\ 

Here L = (Zjc, Ly, Lz) where Lx = 1 I , Ly = I and Lz — —10 . L is the angular momentum; 

\0 -1 0/ \ 1 / \ 0/ 

it is odd under T and even under P. rip = 4= (1,0, 1) is perpendicular to the T-N-P-H plane where all the degeneracies lie in. 

rio = 4j(— 1,0, 1) is along T-N. B is the external magnetic field. (X\, a 2 , 71, y 2 and j3 are positive real parameters that control 
the relative strengths of the different terms in the Hamiltonian. 

The original Hamiltonian Ho[l] produces the three-fold quadratic degeneracy at zero energy, same as the PhC bandstructure 
at r in Fig. 2a. ot\ controls the curvature of the top band; a 2 controls the curvatures of the two degenerate bands concaving 
down. 

Perturbation V\ breaks neither P nor T symmetries. This perturbation lifts the three-fold degeneracy at T while keeping the 
quadratic degeneracy between the lower two bands at zero energy: this reproduces the linear line-degeneracy result in Fig. 2b. 

Perturbation y 2 breaks P. Vi and y 2 together reproduce the two pairs of Weyl points along T-N and T-H in Fig. 2c. 

Perturbation V3 breaks T. When this perturbation is turned on alone, a pair of Weyl points of opposite chiralities emerge along 
the direction of the magnetic field (B), which is consistent with our PhC result in Fig. 2d. 

When all three perturbation terms are non-zero, this model reproduces the same II-IV phase diagrams in Fig. [A] (i.e. when B 
is along T-H). 

ACCIDENTAL LINEAR-DEGENERACIES 

2D or 3D linear dispersions of two bands can also exist when being accidentally degenerate among other bands at the Bril- 
louin zone(BZ) center; but those points do not possess any topological features ||2|-|4j], so they are quite different from the work 
presented here. 
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